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UNDERSTANDING COMPACT OBJECT FORMATION AND NATAL KICKS. IV THE CASE OF IC 10 X-l 
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ABSTRACT 

The extragalactic X-ray binary IC 10 X-l has attracted attention as it is possibly the host of the most massive 
stellar-mass black-hole (BH) known to date. Here we consider all available observational constraints and 
construct its evolutionary history up to the instant just before the formation of the BH. Our analysis accounts for 
the simplest possible history that includes three evolutionary phases: binary orbital dynamics at core collapse, 
common envelope (CE) evolution, and evolution of the BH-helium star binary progenitor of the observed 
system. We derive the complete set of constraints on the progenitor system at various evolutionary stages. 
Specifically: right before the core collapse event, we find the mass of the BH immediate progenitor to be 
> 31 Mq (at 95% of confidence, same hereafter). The magnitude of the natal kick imparted to the BH is 
constrained to be < 130 km/s. More significantly we obtain robust constraints on the physics of the CE event. 
IC 10 X-l provides leverage in distinguishing various treatments of CE evolution: we find that the "enthalpy" 
formalism recently suggested by Ivanova & Chaichenets is the only one that can explain the existence of IC 10 
X-l without the need of invoking unreasonably high CE efficiencies. With this physically motivated formalism, 
we find that the CE efficiency required to explain the system is in the range of ~ 0.6-1. 
Subject headings: binaries: close — stars: evolution — X-rays: binaries — X-rays: individual (IC 10 X-l) 



1. INTRODUCTION 

Over the past few decades, it has become clear that neu- 
tron stars (NS) receive recoil kicks at birth (also known 
as natal kicks) during the core collapse event. This con- 
clu sion is based on proper motion studies of pulsars (see 
e.g. Gunn & Ostriker 1970; Lvne et al. 1982; Lvne & Lorimer 
1994; Brisken et al. 2003; Hobbs et al. 2005; Chatterjee et al. 
2009) and evolutionary studies of NS-hosting binaries 



(see e.g. [ Brandt & Podsiadlowskil 119951; iPfahl et all 



2002; 



2009 



Thorsett et al. 2005; Willems et al. 2006; Martin et al. 
Wong et al. 2010). However, whether black holes (BH) re- 
ceive similar natal kicks during the co r e collapse event is still 
uncertain (see e. g. Brand t et all 1 19951: [Nelemans et all 1 19991; 
Gualandris eTail 120051: [Dhawan et alj I2007t iRepetto et all 



20121) . If BH kicks are ubiquitous, then BH formation must 



be closely associated to that of NS before the formation of the 
event horizon. Otherwise, the formation of BHs through more 
than one physical process will be favored. 

To shed light on questions related to BH formation, we per- 
form detailed evolutionary modeling of the observed BH X- 
ray binaries (XRB), which enables us to derive robust con- 
straints on the mass of the BH immediate progenitor and the 
magnitude of the natal kick imparted to BH during the core 
collapse event. In the previous thr ee papers of this seri es, we 
studied the XRB GRO J 1655-4 (Willems et aD l2005T). XTE 
Jl 118 +480 dFragos et alll2009l) and Cygnus X-l dWong et all 
120121) . For GRO J1655-40, we constrained the mass of the 
BH immediate progenitor to be ~ 5.5-1 1.0 M Q and the mag- 
nitude of its natal kick to be < 210 km/s. In the study of 
Jl 118+480, we found that a natal kick of magnitude ~ 80- 
310 km/s is required to explain the formation of this system, 
and also derived a lower limit of ~ 6.5 M Q on the mass of 
the BH immediate progenitor. Finally, for Cygnus X-l we 



constrained the mass of the BH immediate progenitor and the 
magnitude of its natal kick to be > IS M^ an d < 77 km/s, re- 
spectively. Similarly, Valsecchi et al. (2010) studied the for- 
mation of M33 X-7 using binary modeling as well and derived 
the mass of the BH immediate progenitor and the magnitude 
of the natal kick to be in the range of ~ 15.0-16.1 M and 
~ 10-850 km/s, respectively. In this paper, we study the for- 
mation of the BH in the extragalactic XRB IC 10 X-l. 

IC 10 X-l is one of the three observed BH XRBs that are 
known to host a Wolf-R ayet (WR) star as the mass donor 
dClark & Crowther 2004). Th e othe r two systems are Cygnus 
X-3 (see e. g. ZdziarskiefaI||20T3l) and NGC 300 X-l (see 
e.g. Crowther et al. 2010). The X-ray emission is powered 
by the accretion of stellar wind material onto the BH. At 
the present time, the BH in IC 10 X-l is the most mas- 
sive known stellar-mass BH (^ 23- 34 M , iPrestwich et all 
120071; ISilverman & F ilippenko 2008). Since the supergiant 
progenitor of the observed WR star cannot fit into the 
tight orbit at the present time (orbital period = 34.93 hr, 
Silverman & Filippenko 2008), it is natural to consider the 
system's evolution via a common envelope (CE) evolution 
phase, which involves the BH and the massive progenitor of 
the observed WR star. However, such a binary is likely to 
merge at the end of the CE phase, as the envelopes of massive 
stars are tightly bound (Podsiadlowski et al. 2003). 

iDe Mink et all (2009) suggested an alternative formation 
scenario for IC 10 X-l that does not invoke any CE phase, 
which they called "Case M" evolution. They considered two 
massive stars in a tight orbit, such that tidal interactions al- 
ways kept them spinning rapidly. This led to efficient mixing 
of elements throughout their stellar interiors via rotational ef- 
fects, and hence both stars went through chemical homoge- 
nous evolutions. They stayed compact throughout their main 
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sequence evolution and turned into abnormally massive he- 
lium (He) stars without going through any CE evolution. Al- 
though Case M helps to explain the high masses of the BH 
and the WR star in IC 10 X-l, the short orbital period gives 
rise to a difficulty: the intense mass loss suffered by the WR 
star and its progenitor will widen the orbit. It is very hard to 
explain the current tight orbit without a CE phase. As a result, 
the evolution history of IC 10 X-l has re mained uncer t ain. 

Instead of studying its past evolution, iBulik et ail (120111) 
performed binary modeling to predict the fate of IC 10 X- 
1 . They estimated that the observed WR star will go through 
core collapse in < 0.3 Myr, leading to the formation of a close 
double BH binary with a short coalescence time (~ 3 Gyr). 

The plan of the paper is as follows. In ^2] we review the 
current available observational constraints of IC 10 X-l. A 
general outline of our analysis methodology is presented in 
^3j while detailed discussions of each individual step are in 
$4]-[6] Our derived constraints related to the formation of the 
BH and the past evolution of IC 10 X-l are discussed in $2l 
and we offer our conclusions in the final section. 

2. OBSERVATIONAL CONSTRAINTS FOR IC 10 X-l 

IC 10 X-l is a persistent X-ray sour ce in the local star- 
burst galaxy IC 10. It was discovered by Brandt et al. ( 1997) 
in their X-ray observations of IC 10 made with ROSAT 
iBauer & Brandt! d2004l) derived the 0.5-8.0 keV unabsorbed 
X-ray luminosity of IC 10 X-l to be 1 .50 x 10 38 ergs/s from 
their Chandra observations. The optical counterpart of IC 
10 X-l was ide ntified as the WR star [MAC 921 17A by 
IClark & Crowthe r (2004). Based on the ISchaerer & Maederl 
( 1992) mass-luminosity relationship, they estimated the mass 
of this WR star to b e 35 M w . Using the da ta from Chandra and 
Swift observations. iPrestwich et all (12007b determined the or- 
bital period of IC 10X-1 to be 34.40 ± 0.83 hr and obtained a 
mass function of 7.8 M Q . Using their mass function and con- 
sidering an inclination angle of 90 degrees, they estimated the 
mass of the BH to be 23-34 Mq. Their results indicated that 
IC 10 X-l hosted the most massive stellar mass BH known 
at the present time. The observed orbital period also im- 
plied that the BH is currently accreting mass from the intense 
stellar wind of its WR companion. ISilverman & Filippenko 
(2008) precisely measured the radial velocity amplitude of the 
WR star and refined the orbital period and mass function to 
34.93 ± 0.04 hr and 7.64 ± 1.26M , respectively. For conve- 
nience, our adopted observational constraints are summarized 
in TableQ] 

3. OUTLINE OF ANALYSIS METHODOLOGY 

Among the evolutionary scenarios that could potentially ex- 
plain the formation of IC 10 X-l, we adopt the simplest pos- 
sible evolutionary history that provides a consistent explana- 
tion of all observational constraints. First of all, we assume 
that the BH progenitor and its companion were born at the 
same time. Towards the end of the BH progenitor's life, it 
lost its hydrogen (H) rich envelope because of mass loss via 
a stellar wind or binary interactions. Hence, the BH immedi- 
ate progenitor is a helium (He) rich star. Soon after the birth 
of the BH, the companion star evolved off the main sequence 
and became a supergiant. Eventually, the star overfilled its 
Roche lobe and underwent a phase of dynamically unstable 
mass transfer, which inevitably led to CE evolution. During 
this CE phase, the BH was engulfed into the H-rich envelope 
of its companion. The CE phase ended when the envelope was 
ejected, leaving behind a binary consisting of a BH and a He 



star in a very tight, circular orbit. The intense mass loss via 
stellar wind from the He star continuously drove the binary 
components further away from each other. A very small frac- 
tion of this stellar wind material was accreted onto the BH, 
resulting in X-ray emission from IC 10 X-l that we see today. 

This is the simplest possible evolutionary channel one can 
envision given the current properties of the system. Therefore 
in this paper, we restrict ourselves to the formation of IC 10 
X-l through the above evolutionary channel. In order to de- 
rived constraints on the formation of the BH in IC 10 X-l, we 
track the evolutionary history of the system back to the instant 
just before the core collapse event. Our analysis incorporates 
a number of calculations that can be grouped in three main 
steps. Hereafter, we add the prefix "pre-" and "post-" to the 
name of any event occurred in the evolutionary history of IC 
10 X-l, in order to indicate the instant just before and right 
after that event, respectively. 

In the first analysis step, our goal is to derive the post-CE 
binary and stellar properties of IC 10 X-l. We start by con- 
structing He star models that satisfy the current mass and lu- 
minosity constraints given in Table Q] Using the properties 
of our He star models and different BH masses and post-CE 
binary semi-major axes, we evolve the post-CE binary's or- 
bit forward in time until the current epoch. This calculation 
accounts for tides, wind mass loss, wind accretion onto the 
BH and orbital angular momentum loss via gravitational radi- 
ation. We examine the evolutionary sequence of every binary 
model to find whether at the present time the BH mass and the 
orbital period simultaneously match the observational data, in 
which case we classify that evolutionary sequence as "suc- 
cessful". Furthermore, we consider 1 a of observational un- 
certainties when matching the properties of our models with 
the observational data. The post-CE binary properties can 
then be obtained from our successful evolutionary sequences. 

In the second analysis step, we study the CE event to deter- 
mine whether our derived properties of the post-CE binary can 
be achieved. To do so we employ the standard a-prescription 
(Webbink 1984), examining also several alternative formula- 
tion of the energy budget prescription, to compute the CE ef- 
ficiency, acE- If this parameter is < 1, we conclude that our 
derived properties of the post-CE binary can be explained by 
the current understanding of CE evolution. 

In the final analysis step, we utilize the results from the first 
two steps and perform a Monte Carlo simulation of the orbital 
dynamics involved in the core collapse event. Our goal is to 
derive limits and construct probability distribution functions 
(PDF) for the BH immediate progenitor mass and the poten- 
tial natal kick magnitude imparted to the BH. We start with 
randomized properties of the pre-core collapse (pre-SN) bi- 
nary, which include BH immediate progenitor mass, orbital 
semi-major axis, eccentricity and phase. The core collapse 
event is approximated as occurring instantaneously, with the 
mass ejection and possible asymmetries in the explosion im- 
parting to the BH a potential kick that has random magnitude 
and direction. Using the equations of orbital energy and an- 
gular momentum, we determine whether the binary can sur- 
vive the core collapse event. If it does, we map the pre-SN 
binary properties to the post-core collapse (post-SN) phase 
space. Then, we apply the constraint related to CE evolution 
and retain the data points with olqe < 1 . At the end, our con- 
straints related to the BH formation in IC 10 X-l are derived 
from the retained data points. 

To build the stellar models necessary for our analysis, 
we use a version of the STARS code originally devel- 
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Table 1 

Properties of IC 10X-1 



Parameter 



Notation 



Value 



References 



Distance (kpc) 

Orbital period (hr) 

Inclination angle (deg) 

Mass Function (Mq ) 

Black hole mass (Mq) 

Companion mass (Mq) 

Companion Luminosity (10 6 L@) 

Companion Effective temperature (K) 

0.5-8.0 keV unabsorbed X-ray luminosity (10 38 erg s ) 

Metallicity (Z ) 



d 

Porb 

/Waff) 

M B h 

M 2 

L 2 

T 'ff2 

Lx 

Z 



590 ±35 
34.93 ±0.04 

65-90 

7.64 ±1.26 

~23-34 

17-35 

1.122 x(rf/590 kpc) 2 

85 000 

1.50 x(rf/700 kpc) 2 

0.2 



(1), (3) 
(6) 

(6) 
(6) 
(5) 
(3), (5) 
(3) 
(3) 
(2) 
(4) 



References. — (1) Borissova et al. 2000; (2) Bauer & Brandt 2004; (3) Clark & Crowther 2004; (4) Leroy et al. 2006; (5) Prestwich et al. 2007; (6) Silverman 
& Filippenko 2008. 



oped by Peter Eggleton (Eggleton 1971, 1973; Pols et al. 
1995; Eggleton & Kiseleva-Eggleton 2002; Glebbeek et al. 
2008). The adopted opacity tables co mbine the OPAL opac- 
ities from Iglesias & Rogers (1996), the low temperature 
molecular opacities from Ferguson et al. (2005), the electron- 
conduction opacities from Cassisi et al. (2007) and the Comp- 
ton scattering opacities from Buchler& Yueh (1976). The 
assumed heavy element composition in each stellar model 
is scaled according to th e solar abundances described by 
Anders & Grevesse (1989). Mixing of elements due to con- 
vection is taken into account (Eggleton] 119731) . but the ef- 
fects of rotational mixing and meridional circulation are ex- 
cluded. We notice that these effects can be very impor- 
tant. However, their efficiencies are poorly constrained and 
need to be tuned using the observational data. We choose 
not to add extra free parameters because our models are in- 
tended to be the simplest ones able to satisfy the observa- 
tional constraints. For mass loss during the main sequence 
evolution, we adopt the prescriptions of Vink et al. (2001). If 
the surface H abundance drop s below 0.4, we s witch to the 
WR prescription developed by Nugis & Lamers (2002) with 
the m etallicity scaling law determined by Vink & de Koter 
(2005). When the effective temperature falls below the work- 
ing regime (i.e. < 10 kK) of the mass loss prescriptions de- 
veloped by Vink et al. (2001), we use the mass loss rates de- 
termined by yeJagereFal. (1988) instead. 

4. POST-COMMON ENVELOPE BINARY MODELING 

To model the post-CE phase of IC 1 X- 1 , we start by mod- 
eling the observed WR star. We construct He star models with 
different initial masses at the measured metallicity for IC 10 
(Z = 0.004, see Table [TJ using the stellar evolutionary code 
described in §|3] Then, we evolve the models and retain those 
that at a certain time during their evolution have a luminos- 
ity and mass in agreement with the observations listed in Ta- 
bleQ](see Figure[l]i. Hereafter, we refer to these He star mod- 
els as "successful". At this stage, we use the observational 
constraint on the WR mass rather than effective temperature. 
This is because the effective temperature shown in Table [TJis 
indeed the spectroscopic temperature, which is generally un- 
reliable for WR stars (see Crowther 2007). 

Next, we use the properties of the successful He star mod- 
els to follow the evolution of the post-CE binary orbit and the 
He star's spin until the present time. Our goal is to constrain 
the post-CE binary properties. For each successful He star 
model, we vary the post-CE orbital semi-major axis in steps 
of 0. 1 Rq and consider different BH masses according to the 
measured mass function (/(Mbh)) and inclination angle (/) 
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Figure 1. Luminosity as a function of mass for He star models with metal- 
licity Z = 0.004 and initial masses ranging between 32.2 and 41.6 Mq. The 
grey region represents the observational constraints shown in Table[T] 

listed in Table [1] Specifically, we take the mass of each suc- 
cessful He star model at the current epoch and compute the 
maximum and minimum BH mass at present (MBH,max and 
MBH,min, respectively) using the measured /{Mbh) and i with 
lcr of observational uncertainties. Then, we define 



M, 



BHjnean '- 



am r 



Mbh. max + MBH,min 



Mbh, max- Mbh. i 



(1) 

(2) 



and vary the BH mass at present from MBH.min to Mbh. max, 
in steps of 0.2 • AM BH . The amount of mas s accreted onto 
the BH since its birth is negligible, as the Bo ndi & Hovld 
(1944) accretion of the stellar-wind material leads to a very 
small capturing fraction. Hence, the post-CE BH mass is only 
slightly different from the one at present. For the post-CE 
configuration of the binary orbit and the He star's spin, we 
assume that CE evolution circularizes and synchronizes the 
binary orbit. We note that our assumption that the post-CE 
He star's spin is synchronized with the orbital frequency is 
completely arbitrary. However, this arbitrary choice only af- 
fects the strength of the tide exerted on the He star by the 
BH. S ince this tide is usually weak (e.g. see iDetmers et alJ 
2008), our arbitrary assumption of the post-CE He star's spin 
does not have any significant influence on our analysis. Under 
these assumptions, we follow the evolution of the binary or- 
bital semi-major axis and eccentricity, and the He star's spin, 
accounting for tides exerted on the He star by the BH, stellar 
wind mass loss, orbital angular momentum loss due to gravi- 
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Figure 2. The variation of our predicted X-ray luminosity against the BH 
mass at present for every successful evolutionary sequences. We consider two 
extreme values for the spin of the BH (see Equation|5]: non-spinning (black 
dots) and maximally spinning (orange dots). The grey region represents the 
observed X-ray luminosity presented in Table [T] 



tational radiation, and accretion of the He star's wind material 
onto the BH. In order to determine how much wind mate- 
rial is accreted onto the BH, we compute the accretion rate 

Ma, 



Specifically, we adopt the Bondi & Hoyle (1944) ac- 
cretion model and follo w the formalism in § 4.1 and 4.2 of 
iBelczvnski et al.l d2008l) to obtain M acc . 

The relevant ordinary differential equations (ODEs) gov- 
erning the orbital evolution are integrated forward in time. 
To do this, we use the orbital evolution code described in the 
supplementary information of Valsecchi et al. (2010) with the 
following two modifications. 

1 . For the second- order tidal coe fficient E%, we take a stel- 
lar model from IClaretl (|2005) with an initial mass of 
25.2 M & and metallicity of 0.004 (both in agreement 
with the observations shown in Table [TJ, and derive 



log(£ 2 ) =-5.49491 
-1.94284? 



51.277 



1.99707? 



2.41139 



(3) 



where tuis is the star's evolutionary time expressed in 
units of the main sequence lifetime. As we are dealing 
with a He star, tyis is taken to be the time the star has 
spent burning He at its center. 

2. For the effect of wind mass loss on the spin of the He 
star, we approximate that the wind carries away its an- 
gular momentum from a thin shell at the star's surface. 
Hence, we set 



J spin = ~M 2 R\u, 



(4) 



where J sp j„ is the spin angular momentum, R 2 is the 
stellar radius and u is the spin angular frequency of the 
He star . 

For each combination of post-CE binary component masses 
and orbital semi-major axis, the integration of the relevant 
ODEs proceeds forward in time only if the He star is spin- 
ning slower than the break-up frequency il c « ^JGMJR? and 
is not filling its Roche lobe. 

At the end of the calculation, we retain the evolutionary se- 
quences that at the present time match the observed orbital 
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Figure 3. Post-CE binary component masses and semi-major axis of IC 10 
X-l given by all successful evolutionary sequences. 



period. Hereafter, we refer to them as "successful" evolution- 
ary sequences. The X-ray luminosity at the present time of 
each evolutionary sequence is computed by following the for- 
malism in § 9.1 of Belczvnski et al. (2008). Using the relation 
between the radius of the accretor and the unknown BH spin, 
the X-ray luminosity is given by 



Lx = Vboie 



M BH c 2 



(5) 



where c is the speed of light, and rjb„i and e are adopted to 
be 0.8 and 0.5, respectively. The variable r; sco is the radius 
of inner most stable circular orbit around the BH expressed 
in units of GMbh/c 1 . It equals 6 for a non-spinning BH 
and 1 for a maximally spinning BH. Figure [2] shows that for 
every successful evolutionary sequence the observed X-ray 
luminosity always falls in between the upper (r,-, co =1) and 
lower (fisc = 6) limits of our predicted Lx at the present time. 
Hence, we are not formally imposing the observed Lx as a 
constraint, but our predicted Lx at the present time of all suc- 
cessful evolutionary sequences are naturally consistent with 
observed values. 

The post-CE binary component masses and orbital semi- 
major axis of IC 10 X-l given by all successful sequences are 
illustrated in Figure [3] The figure shows that the post-CE bi- 
nary consists of a 24.8-38.7 M Q BH and a 32.2-41. 6 M He 
star and has an orbital separation of 17.1-22.3 Rq. The equiv- 
alent orbital period is 25.3-35.0 hr. Due to the intense stellar 
wind suffered by the He star, the orbital separation increases 
continuously. Such a binary evolves to become the observed 
XRB IC 10 X-l in < 0.36 Myr. 

5. COMMON ENVELOPE EVOLUTION FORMALISM 

Prior to the Monte Carlo simulation of the orbital dynam- 
ics involved in the core collapse event, we study whether our 
derived post-CE binary properties of IC 10 X-l is achiev- 
able based on the current understanding of CE evolution. We 
construct stellar models of pre-CE He star progenitor candi- 
dates, which have core masses covering the entire range of 
post-CE He star's mass given by the successful evolutionary 
sequences. These models are built based on the stellar evolu- 
tion of isolated stars. We note that the BH progenitor could 
have potentially transferred mass to its companion before the 
core collapse event. However, we assume the He star progen- 
itor will quickly adjust itself after the end of mass transfer. 
Hence, its structure right before the CE phase is not expected 
to be significantly different from an isolated star of the same 
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Figure 4. Variation of maximum AE or b against the mass of thepost-CE 
He star (Mwr)- The maximum A£ or j is computed with Equation (7), under 
the assumption that the pre-CE semi-major axis (A„ re cg) approaches infinity. 
The post-CE parameters are given by successful evolutionary sequences. 

mass. Furthermore, in order to initiate CE evolution, the im- 
mediate progenitor of the He star needs to fill its Roche lobe 
and start unstable mass transfer before its radius reaches its 
maximum value at the instant t Rmax . If this is not the case, the 
star will shrink rapidly afterward due to intense stellar wind 
mass loss. Meanwhile, most of this wind material will leave 
the binary system, resulting in the widening of the binary or- 
bit. Because of the increase in orbital separation, the Roche 
lobe of the star expands. Since the star is shrinking and its 
Roche lobe is expanding, it cannot overfill its Roche lobe and 
initiate CE evolution after t Rmax . Thus, we only consider the 
properties of our He star progenitor models before t Rmax when 
studying the CE event. 

To examine whether CE evolution can explain our derived 
post-CE properties of IC 10 X-l, we use the energy formal- 
ism ( Webbink 1984) and compute the corresponding common 
envelope efficiency acE, 

acE ■ AE or b = Ebind , (6) 

where AE or i, is the change in orbital energy and Eund is the 
energy required for dispersing the envelope to infinity. Here, 
AE or j, is simply 



GM BH M WRpw GM BH M WR 

L±E or b = — 1- ■ 



2A 



preCE 



2A 



(7) 



postCE 



where Mwr P w and Mwr are the masses of the pre-CE He star 
progenitor and the post-CE He star, and A preCE and A pmtCE 
are the pre- and post-CE orbital semi-major axis, respectively. 
Since the outcome of CE phase is given by the successful 
evolutionary sequences, the second term in Equation (|7]l is 
predetermined. Hence, it is obvious that AE or b reaches maxi- 
mum when A pre cE approaches oo. Depending on the post-CE 
parameters given by each successful evolutionary sequence, 
maximum AE or b ranges from 7.30 x 10 49 to 1.58 x 10 5() erg 
(see Figure|4]i. 

When computing Ebmd of our He immediate progenitor 
models, we set the core-envelope boundary at the base of the 
convective thick hydrogen burning shell. Since this is gen- 
erally located close to where the hydrogen abundance (Xh) 



is 10%, we si mply defin e the core boundary at Xh = 0.1. 
iDewi&Taunsl feOOO) and lTauris & Dewil (12001b also consid- 
ered the same definition of core boundary in their studies of 
CE evolution. 

Based on our adopted energy formalism shown as Equa- 
tion |6j, olqe < 1 means the loss in the orbital energy is suffi- 
cient to unbind to envelope of our He star immediate progeni- 
tor model. In order words, this means CE evolution is capable 
of producing our derived post-CE properties of IC 10 X-l. On 
the other hand, chqe > 1 means that there is not enough energy 
to unbind the envelope and the pre-CE binary will eventually 
merge instead. We exa mine four different treatments o f CE 
evo lution: t he original Webbink (1984) prescription (§ 15. U . 
the Webbink ( 1984 ) prescription with enhanced convective 
overshooting ( §I5.21 i. allow ing h yper-critical accretion onto the 
BH during CE evolution ( §15.31 ). a nd th e "enthalpy" formalism 
(llvanova & Chai chenets 2011) ( §54i i. We find that only the 
"enthalpy" formalism is capable of explaining the post-CE bi- 
nary properties given by all successful evolutionary sequences 
with acE < 1. In the following subsections, we will discuss 
each of the aforementioned treatments of CE evolution. 



5.1. Orisinal Webbinli A1984) Prescription 
First, we use the standard definition of Ebind, which can be 



written as 



/•surface 

Ebind =~ ($(m) + e(m)) dm. 

J core boundary 



(8) 



Here $(m) = -Gm/r is the gravitational potential and e(m) 
is the specific internal energy. In Figure we illustrate the 
typical variation of Ebind throughout the evolution of a He star 
progenitor. After the termination of main-sequence at t tms , 
Ebind increases as the He star progenitor shrinks. It reaches 
a maximum shortly after the beginning of convective thick H 
shell burning at tg S b- Then, it starts to decrease because of the 
envelope expansion. Figure [5] also shows that throughout its 
evolution prior to t Rmax , Ebind of a He star progenitor is always 
the smallest at t Rmax . 

Using the above definition of Ebmd and Equation (O, we 
search for He star progenitors models that will give us acE < 
1 at a certain time during its evolution before t Rmax . To achieve 
this goal, we create a grid of stellar models at the observed 
metallicity of IC 10 X-l by varying the initial mass from 65 to 
105 Mq in steps of 0.1 M@. This grid covers the entire initial 
mass range of massive stars, which can have core masses fall 
within the mass range of post-CE He stars constrained by all 
successful evolutionary sequences (i.e. 32.2-41.6 Mq). 

For each model in our grid, we match its core mass at any 
time before t Rmax with the post-CE He star masses of the suc- 
cessful evolutionary sequences. Then, we compute AE or b 
with Equation (J7]) using the relevant post-CE binary proper- 
ties given by those matched sequences. Using this AE or b and 
the Ebind defined by Equation ([8]), we can obtain acE with 
Equation ©. We find that for all models oiqe > 1, as Ebmd is 
always several times larger than AE or b- To illustrate this, we 
select six representative models of He star progenitors from 
our grid and list them in Table|2]as Model I-VI. Here, we con- 
sider the CE evolution involving these models at t Rmax , which 
is the time when these models have the lowest Ebind- The mass 
range of the cores in these models at t Rmax is the same as that 
of the post-CE He stars given by all successful evolutionary 
sequences. In column 8, we list the maximum AE or b during 
the CE evolution involving these models, under the assump- 
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Table 2 
Properties of Selected He Star Progenitor Models 



Model 


Initial Mass (Mq) 


a ov a 


MwRpro (M Q f 


M core (M @ Y 


RwRpro (Rq) 


E bmd (10 50 erg) f 


max. AE orb (10 50 erg/ 


I 


80.4 


0.12 


65.9 


32.2 


1406 


4.05 


1.14 


II 


84.1 


0.12 


67.3 


34.1 


1337 


4.33 


1.23 


III 


88.6 


0.12 


69.9 


36.1 


1268 


4.72 


1.32 


IV 


92.6 


0.12 


72.9 


38.2 


1214 


5.13 


1.42 


V 


95.9 


0.12 


75.6 


40.0 


1175 


5.48 


1.50 


VI 


99.9 


0.12 


78.2 


41.6 


1136 


5.85 


1.58 


VII 


75.0 


0.20 


60.0 


32.2 


1429 


3.84 


1.14 


VIII 


78.5 


0.20 


62.4 


34.1 


1343 


4.14 


1.23 


rx 


83.0 


0.20 


66.3 


36.1 


1262 


4.63 


1.32 


X 


86.8 


0.20 


68.2 


38.2 


1205 


4.85 


1.42 


XI 


89.9 


0.20 


70.9 


40.0 


1171 


5.16 


1.50 


XII 


93.6 


0.20 


73.0 


41.6 


1140 


5.46 


1.58 


XIII 


67.0 


0.30 


53.8 


32.2 


1403 


3.49 


1.14 


XIV 


70.9 


0.30 


56.9 


34.1 


1317 


3.82 


1.23 


XV 


74.8 


0.30 


59.7 


36.1 


1239 


4.13 


1.32 


XVI 


78.6 


0.30 


62.4 


38.2 


1188 


4.40 


1.42 


XVII 


81.8 


0.30 


65.1 


40.0 


1149 


4.69 


1.50 


XVIII 


84.8 


0.30 


67.0 


41.6 


1118 


4.92 


1.58 



1x10 3 



9 x 1 3 



3x1 cr 



7x1 3 



6xi cr 



5x10= 



Model II 
Model IX 
Model XVIII 




Note. — The models with a ov = 0. 12 have normal strength of convective core overshooting. The parameters listed in column 4-7 are the properties of the He 
star progenitor models when their radii reach their maximum values at the instant tn,„ ax . 
" convective overshooting parameter 
* mass of the He star immediate progenitor 

c core mass of the He star immediate progenitor (same as the mass of the descendent He star right after CE evolution) 
d radius of the He star immediate progenitor 

e energy required to disperse envelope of the He star immediate progenitor 
f maximum change in orbital energy during CE evolution involving the He star progenitor at tsmu ( see Figure[4j 

star progenitor models with similar core masses but less mas- 
sive envelopes. We cannot simply move our core boundary 
to a higher Xh- This is because doing so will include part of 
the convective thick H burning shell as core once its Xh drops 
below the newly chosen value, which contradicts with our def- 
inition of core-envelope bifurcation point. Instead, we choose 
to increase the convective core overshooting {a ov ) parameter 
from a normal value of 0.12 to 0.2 and 0.3. In this trial, we 
again compute uce using Equation (O with Ebmd defining as 
Equation ©. 

As before, we make grids of He star progenitor models at 
the observed metallicity of IC 10 X-l with a ov = 0.2 and 0.3, 
by varying the initial mass between 60-95 and 55-88 Mq in 
steps of 0.1 Mq, respectively. Each grid covers the entire ini- 
tial mass range of massive stars, which can have core masses 
fall within the mass range of the post-CE He stars constrained 
by all successful evolutionary sequences (i.e. 32.2-41.6 Mq). 

Aft er co mputing olqe in the same way as our previous trial 
(see 35. lb for all He star progenitor models in our grids of 
a ov = 0.2 or 0.3, we do not find any cases of acE < 1. To 
help explaining this result, we select a group of representa- 
tive models from each grid and list them as Model VII-XII 
(a ov = 0.2) and Model XIII-XVIII (a OT , = 0.3) in Table |2| 
Here, we consider the CE evolution involving these mod- 
els at tRmax, which is the time when these models have the 
lowest Ebind (see Figure |5). For each group of models, their 
core masses cover the entire range of post-CE He star's mass 
given by the successful evolutionary sequences. The maxi- 
mum AE or b during the CE evolution involving these models 
(see Figure 2|l are listed in column 8. When comparing the 
Ebind of these models to that of the models with the normal 
strength of convective overshooting (i.e. Model I- VI), we 
find that for the models with the same core mass Eb m d only 
decreases slightly with the increasing a ov . Considering the 
difference between Ebmd and maximum AE or b of these mod- 
els, that amount of decrease in Ebmd is not enough to make 
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log(R«,Rpro/ R o) 

Figure 5. The variation of Ei,j„j against log(RwRp m ) throughout the evo- 
lution of three different He star progenitor models (Models II, IX, XVIII as 
listed in Table [1), which have similar initial masses but different strength of 
convective overshooting. Different events in stellar evolution are illustrated 
on each curve: the termination of main-sequence (triangles), the beginning 
of thick H shell burning (circles) and the time of their stellar radii reaching 
maximum(squares). Throughout the evolution of these models, their £(,,,„/ 
have similar behavior and reach the minimum when the stellar radii of these 
models reach maximum. 

tion that A pre cE approaches oo (see Figure |4). It is obvious 
that for these models Ebmd is at least 3.5 times larger than the 
corresponding maximum AE or i,. Thus, olce is always > 1. 

5.2. Enhanced Convective Overshooting 

A possible reason for the negative result in the previous trial 
could be that the envelopes of our He star progenitor models 
are too massive. Hence, we might obtain olqe < 1 using He 
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a CE < 1. 

5.3. Hyper-Critical Accretion 

It has been suggested that a compact object might accrete 
a significant amount of mass after being engulfed into the 
envelope of its companion, due to hyper-critical a ccretion 
JBlondidl 19861: lChevaheri[l98l fl993t lBrownlfl995h . Since 
only part of the companion's envelope will be dispersed to in- 
finity, the Ei,i n d in the energy formalism shown as Equation © 
needs to be adjusted accordingly. Hence, we write the equa- 
tion of energy balance as 

acE ■ AE or b = AEbind = fej • Eb 



iiul • 



where 



fej 



M„ 



-AM R 



M en 



(9) 



(10) 



is the fraction of envelope mass ejected to infinity. Here, M em , 
is the mass of the envelope and AjMbh is the amount of mass 
acc reted onto the BH . Based on Equation (O and ( fTUl ), we fol- 
low |BelczynskreFal] d2002l) and derive the rates of change in 
the BH mass and the binary semi-major axis during the phase 
of hyper-critical accretion. The detailed derivation of these 
rates with respect to the mass of the He star progenitor (M com ) 
during the accretion phase can be found i n the Appendix. 
Our rates are different from those derived by Belczvns ki et al.l 
(2002), since we consider a fraction instead of the total Ebmd 
when balancing the energy budget of envelope ejection. In 
other words, we use Equation (0 instead of Equation © to 
incorporate oiqe- 

Using the mod els o f He star progenitors in our constructed 
grids (see j j5.1| and 15. 2\ and considering their properties 
throughout thei r evol ution prior to tR max , we numerically inte- 
grate Equation dAlOt and ( IA1 It from M com = M core to MwR pm 
with different occe between and 1 . We find that the BH could 
have accreted ~ 10 M during CE phase. However, this can 
happen only if the He star progenitor at the onset of the hyper- 
critical accretion phase (i.e. when the BH touches its surface) 
has a radius at least 1.3 times larger than its radius just be- 
fore CE pha se (i.e . when it fills its Roche lobe). According 
to Equation ( IA16I ), there does not exist a pre-CE binary semi- 
major axis leading to such condition. Therefore, we conclude 
that the hyper-critical accretion formalism cannot explain the 
post-CE binary properties given by the successful evolution- 
ary sequences. 

5.4. "Enthalpy" Formalism 

llvanova & Chaichenetsl d2011l) argued that enthalpy should 
be considered when calculating the binding energy of the en- 
velope and introduced 



pi 

^bind ' 



surface 



core boundary 



<1>(ot) + e(m) + 



P(m) 

p(m) 



dm. (11) 



Using Model III listed in Table [2] as an example, we illus- 
trate the behavior of its envelope binding energy computed 
by Equation (0 and ( fTTT ) on Figure [6] It is clear that E' hind is 
always smaller than Ebmd at any time during the evolution of 
Model III. This is because the term P(m)/p(m) is always posi- 
tive everywhere within a star. Including it as an energy source 
lowers the binding energy of the envelope, and hence the CE 
efficiency. Figure [6] also shows that E' bind at tR max is smaller 
than the corresponding maximum AE or f, during the CE phase 
(see Figure |4] and Table |2). According to the Equation ©, 
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Figure 6. The variation of Eti„d and EL. against log(RwRp m ) throughout 
the evolution of Model III. Different events in stellar evolution are illustrated 
on each curve: the termination of main-sequence (triangles), the beginning 
of convective thick H shell burning (circles) and the time of its stellar radius 
reaching maximum (squares). The horizontal line indicates the maximum 
A £ g ,.fc d uring CE evolution involving this He star progenitor model (see text 
in j|5.U . According to Equation l(6}, this model will give us «c£ < 1 at f/;,„ av 
using E' bind as the description of envelope binding energy. 

acE could be < 1 if AE or b during the CE evolution involving 
this progenitor model is close to that maximum value. 

To examine whether the "enthalpy" formalism can explain 
our derived post-CE properties of IC 10 X-l, we consider the 
models of He star progenitors in our constructed grid with the 
norm al strength of convective overshooting {a ov = 0.12, see 
35. lb . Throughout the evolution of each model, we match 
its core mass with the masses of the post-CE He stars given 
by the successful evolutionary sequences. Then we compute 
acE with Equation (0 and ( [TTI ). using the post-CE binary 
properties of the matched sequences. It is not a surprise that 
we find many cases of olqe < 1 ■ 

To summarize, the post-CE binary properties given by 
the successful evolutionary sequences can be explained by 
our understa nding of C E evolution, using the standard acE 
prescription (Webbink 1984) with the "enthalpy" formalism 
(Ivanova & Chaichenets 201 1). IC 10 X-l gives a very strong 
constraint on the physics of CE evolution. It is the very first 
time that a system appears to be capable of distinguishing 
among treatments of CE evolution. The "enthalpy" formal- 
ism appears to be the appropriate description on the energy 
budget of envelope ejection during CE evolution, at least for 
the massive stars involved in the formation of IC 10 X-l. 

6. ORBITAL DYNAMICS AT CORE COLLAPSE 

Using the progenitor properties obtained from the previous 
steps (see $4] and [5]) as constraints, we perform Monte Carlo 
simulations on the binary orbital dynamics involved in the 
core collapse event. Our goal is to derive constraints on the 
properties of the BH immediate progenitor and the magnitude 
of the natal kick imparted to the BH. 

Just before the core collapse event, the binary consists of 
the BH immediate progenitor and its companion star, with an 
orbital semi-major axis A pre sN and eccentricity e pre sN. Note 
that since it is not necessary that the pre-SN progenitor has 
experienced any mass-transfer phase, we cannot assume that 
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the pre-SN orbit is circular. Instead we consider the full range 
of possibilities with eccentric pre-SN orbits. The masses of 
the BH immediate progenitor and its companion are M B u pro 
and Mi, respectively. During the core collapse event, the mass 
loss from the BH immediate progenitor and the potential na- 
tal kick imparted to the BH alter the binary orbital properties. 
Hence, the post-SN orbital semi-major axis and eccentricity 
become A pos ,sn and e pos ,sN, respectively. As we assume the 
companion star is not affected by the instantaneous core col- 
lapse event, its properties remain unchanged. 

Using the equation of binary orbital energy and angular mo- 
mentum, the pre- and pqst-SN binary properties are related as 
(lHillslll983tlWong et al]l2012h : 



v? + vi Hpm + 2v k v BHpm cos e k 

= G(M BH +M 2 ) 



G(M BH +M 2 )A postSN ( 1 - e 2 posrSN ) 



IpostSN 



(12) 



= r 1 \V k sin 2 9 k cos 2 4> k + [sin ip(V BH pw + V k cos k ) 



- V k cos ip sin 9 k sin </> J 



(13) 



Here, V k is the magnitude of the natal kick imparted to the 
BH, while 9 k and <\> k describe its direction in the frame of the 
BH immediate progenitor. Specifically, 9 k is the polar angle 
of the natal kick with respect to the relative orbital velocity 
of the BH immediate progenito r, and <j) k is the c orresponding 
azimuthal angle (see Figure 1 in Kalogera 2000, for a graphic 
representation). The variable r, which is the separation be- 
tween the BH immediate progenitor and its companion at the 
moment of core collapse, can be expressed as 



r - A pre SN{ 1 - e pre SN COS £ pre SN) , 



(14) 



where the pre-SN eccentric anomaly £ pre sN is related to the 
pre-SN mean anomaly M. pre sN as 

M-preSN = £preSN-epreSNSin£preSN- (15) 

The relative orbital speed V B H pn > of the BH immediate pro- 
genitor in the pre-SN binary can be written as 



V B Hpro = 



G(M BHpl 



,+M 2 )(-- 



V) 

re.SN / 



IpreSN 



1/2 



(16) 



Finally, the angle ip is the polar angle of position vector of the 
BH immediate progenitor with respect to its relative orbital 
velocity in its companion's frame, which is related to the pre- 
SN orbital parameters as 



sin^ = 



1 preSN 



' ' e preSN> 



r(2A preSN -r) 



1/2 



(17) 



We start our calculation at the instant just before the core 
collapse event. The properties of the companion star are taken 
from a stellar model in a grid constructed by varying the ini- 
tial mass between 65-105 M Q , insteps of 0.1 M Q . It is indeed 
the same grid of models used in finding the correct formalism 
describing CE evolution, which is capable of explaining the 
post-CE binary prop erti es gi ven by all successful evolutionary 
sequences (see jj5.1l and l5~4l i. To obtain the properties of these 
companion star models at the core collapse event of the BH 
progenitor, we need to approximate when the BH is formed 



in the evolutionary time frame of its companion (t B n). Un- 
der the assumption that the BH progenitor and its companion 
are born at the same time, t B u simply equals to the lifetime of 
the BH progenitor. We adopt that to be the lifetime of a stellar 
model with an initial mass of 150 M Q , which is approximately 
2.90 Myr. Other pre-SN binary properties, namely M B u pro , 
A P reSN, epreSN, -MpreSN, are drawn randomly from uniform dis- 
tributions. Supplemented with a natal kick magnitude V k and 
direction angles (9 k , <p k ) drawn randomly from uniform and 
isotropic distributions, we can obtain A pos ,sn and e pos ,sN from 
the pre-SN binary properties using Equation (Q~2) and (1T3V In 
this calculation, M B h is taken directly from a successful evo- 
lutionary sequence, because the BH in our adopted formation 
scenario of IC 10 X-l has only accreted negligible amount of 
mass since its birth. 

For each combination of companion star model and suc- 
cessful evolutionary sequence, we perform 2000 Monte Carlo 
trials. We only retain the data points that satisfy all of the 
following constraints and classify those data points as "suc- 
cessful". 

1. M B H pm is taken to be less than 60 M . This arbitrary 
upper limit is based on our understanding of the stellar 
evolution of massive stars. We will discuss the impact 
of this limit o n our derived constraints related to the BH 
formation in §17.11 



2. The binary must survive through the core collapse 
event. This means A pos ,sn and e pos ,sN obtained from 
Equation ( [T2] | and (TT~3T > needs to have realistic values: 



IpostSN 



> and < e postSN < 1. 



3. Since the core collapse event of the BH immediate 
progenitor is instantaneous, the separation between the 
pre-SN binary components is the same as that of the 
post-SN binary. This gives a constraint as 

ApreSNi 1 _ epreSN COS EpreSN) 

= A pos tSN( 1 _ SpostSN COS EpostSN) , (18) 

which has to be satisfied with a realistic post-SN eccen- 
tric anomaly: |cos(£ ' PO sisn)\ < 1- 

4. Both pre-SN binary components cannot spin faster than 
the breakup angular velocity Vt c k, ^GM/R 3 . Here, 
we assume their spins are pseudo-synchronized to the 
pre-SN orbital angular velocity. Since the BH imme- 
diate progenitor is expected to be He rich due to po- 
tential binary interaction or intense mass loss via stellar 
wind, we approximate it s radius using Equation (3) in 
iFrver & Kalogera! (11997b . 

5. According to our adopted formation scenario of IC 10 
X-l, both components of the pre-SN binary need to fit 
within their Roche lobes at periapsis. 

6. The BH companion in the post-SN binary must later 
fill its Roche lobe at periapsis before its radius reaches 
its maximum value (see §0, which leads to CE evolu- 
tion. As massive stars evolve at the same time scale, 
the time difference between the formation of the BH 
and the onset of CE event is small. Hence, we assume 
the binary semi-major axis and eccentricity remain un- 
changed within this period of time. The outcome of CE 
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phase is indeed constrained by the post-CE binary prop- 
erties of the corresponding successful evolutionary se- 
quence. When the BH companion fills its Roche lobe at 
periapsis, its core mass needs to match the correspond- 
ing mass of the post-CE He star, with a tolerance of 
1 Mq, This tolerance value is chosen according to the 
initial mass resolution in our grid of post-CE He star 
models, whose properties are used in constructing the 
successful evolutionary sequences (see $4}. Further- 
more, the common envelope efficiency olqe determined 
by the standard a prescripti on (W ebbink 1984) with the 
"enthalpy" formalism (see j j5. A\ must be < 1 . 

7. RESULTS 

The elements presented in the previous sections can now 
be combined to establish a complete picture of how we track 
the evolution of IC 10 X-l backwards in time and derive con- 
straints related to the BH formation in this system. We first 
use the modeling of binary evolution and observational con- 
straints to determine the post-CE binary properties. Specifi- 
cally, our "successful" evolutionary sequences at the present 
time simultaneously match the measured component masses, 
He star luminosity and binary orbital period of IC 10 X-l. 
Then, we search for the correct formalism and treatment of 
CE evolution leading to the formation of IC 10 X-l. We find 
that the standard a prescription (Webbink 1984) with the "en- 
thalpy" formalism (Ivanova& Chaichenets 2011) is capable 
of explaining the post-CE binary properties given by our suc- 
cessful evolutionary sequences. Last, we use our findings in 
the two previous steps as part of the constraints applied on 
a Monte Carlo simulation of the binary orbital dynamics in- 
volved in the core collapse event. Each data point in this sim- 
ulation contains seven random parameters: the BH immedi- 
ate progenitor mass, the pre-SN orbital mean anomaly, semi- 
major axis and eccentricity, the magnitude of the natal kick 
velocity imparted to the BH and two angles describing the 
kick direction. These random parameters are drawn from uni- 
form or isotropic prior distributions. If a data point satisfies 
all constraints mentioned in §|6] we classify it as a "success- 
ful" data point. Our results and derived constraints (at 95.4% 
of confidence) presented in what follows are obtained from 
all successful data points, as well as the PDFs illustrated in 
Figure [7] 

7.1. Core Collapse Constraints 

Just before the core collapse event, we find the BH immedi- 
ate progenitor mass (MgHpro) to be 45. 6^40 M Q and its main- 
sequence companion mass (M2) to be 82.4^- \ Mq. The or- 
bital separation (r) between them is SllO^j^J Rq. During 
the core collapse event, the BH immediate progenitor loses 
< 49.7% of its mass, which is < 20.3% of the total mass in the 
pre-SN binary. Possible asymmetries developed in the core 
collapse event can lead to a natal kick (Vj) of < 128 km/s im- 
parted to the BH. The PDFs of V* and MgHpro are presented in 
Figure [7] 

We notice that the MgHpro PDF contains a plateau between 
38 and 60 Mq, which indicates there is no clear upper limit on 
MsHpro beside the expected maximum value based on the stel- 
lar evolutionary theory. Fortunately, our derived upper limit 
on Vk depends very weakly on our adopted value of maxi- 
mum MsHpm- Figure [8] illustrate the two dimensional joint 
Vu-Mbhpw confidence levels. For MsHpw above 38 Mq, the 
boundaries of the confidence levels are almost perpendicular 



to the Vk axis. This means the Vk PDF remains roughly the 
same even if we omit the data points in a certain mass range 
in this regime. Hence, our derived upper limit on Vk will not 
change significantly if we adopt a lower MsHpm upper limit. 

7.2. Common Envelope Evolution Constraints 

Soon after the formation of the BH, the companion star 
evolves off the main-sequence and expands rapidly. When 
it fills its Roche lobe at the periastron, its mass and radius are 
77.6^| Mq and 1009^ R Q , respectively. Meanwhile, the 
orbital separation at periastron {A per i^ pre cE) is 21 86^0° Rq. 
Then, the binary undergoes dynamically unstable mass trans- 
fer, which leads to CE evolution. We find the common enve- 
lope efficiency (cxce) to be > 0.579, which is similar to the 
values determined from studi es of white dwarf b inaries (see 
Nele mans & Tout! 120051: iDe Marco efail 1201 U iDavis et alj 
120121) . However, a direct comparison is not appropriate, as 
both the properties of the examined systems and the assump- 
tions in the analyses are significantly different. 

At the end of CE event, the binary consists of a 24.8- 
38.7 Mq BH and a 32.2-41.6 M Q He star. They are orbiting 
each other in a very tight and circular orbital, which has a ra- 
dius of 17.1-22.3 Rq. Unlike other limits presented in this 
section, the limits on the post-CE binary properties enclose 
the full range of the derived constraints, which are obtained 
from the post-CE binary modeling discussed in §|4] 

8. CONCLUSION 

In this analysis, we track the evolution of IC 10 X-l back- 
wards in time up to the instant just before the core-collapse 
event and study the formation of the BH in this system. This 
covers the following evolutionary phases: binary orbital dy- 
namics at core collapse, CE evolution, and evolution of the 
BHDHe star binary progenitor of the observed system. We 
first focus on the latter and use the modeling of binary evolu- 
tion to construct successful evolutionary sequences, and use 
them to determine the post-CE binary properties. These se- 
quences are referred as successful because their properties 
at present matches these observational constraints of IC 10 
X-l: binary orbital period, component masses and luminos- 
ity of the WR star. Our predicted X-ray luminosity at the 
present time, resulting from the stellar-wind accretion onto 
the BH, is consistent with the observed values. We then an- 
alyze the evolution through the necessary CE phase. We em- 
ploy different CE treatments, as the standard treatment leads 
to unphysical results. We find that only the "enthalpy" for- 
malism (Ivanova & Chaichenets 201 1) along with an energy- 
based CE efficiency (Webbink 1984) can explain physically 
the post-CE binary properties of the IC 10 X-l progenitor. 
Finally, we perform a Monte Carlo simulation on the or- 
bital dynamics involved in the core collapse event. Each data 
point contains seven free parameters drawn from uniform and 
isotropic distributions, which describe the properties of the 
pre-SN binary and the natal kick imparted to the BH. Besides 
the constraints related to the core collapse event, we also use 
what we learned about the CE event involved in the formation 
of IC 10 X-l as additional constraints to reject data points. If 
a data point satisfies all the constraints mentioned in §|6] such 
as the survival of the binary through the core collapse event 
and the common envelope efficiency oice being < 1, we clas- 
sify it as a successful data point. Our constraints (at 95.4% of 
confidence) related to the BH formation and the CE event oc- 
curred in the past of IC 10 X-l are derived from all successful 
data points. 
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Figure 7. Probability distribution functions (PDF) of different variables (from left to right): (top row) the natal kick magnitude (Vi) imparted to the BH, the BH 
immediate progenitor mass (MgHpro), the pre-CE orbital separation at periastron (A™.; t preCE)< (bottom row) the pre-CE He star progenitor's mass (MwRpm) and 
radius (RwRpm), and the common envelope efficiency (cxce)- 
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Figure 8. Two dimensional joint Vk-^BHpm confidence levels: 68.3% (red), 
95.4% (yellow), and 99.7% (blue). 



We find that in the pre-SN binary, the BH immediate 
(He rich) progenitor has a mass of 45 .6^40 M®, with a 
82.4+7 j Mq main-sequence companion. During the core col- 
lapse event, the BH immediate progenitor loses < 49.7% of 
its mass, which is < 20.3% of the total mass in the pre-SN bi- 
nary, and leaves behind a 24.8-38.7 M Q (full range) BH. The 
asymmetries developed in the core collapse event can lead to 
a natal kick of < 128 km/s imparted to the BH. For the BH 
immediate progenitor mass, we do not find any absolute up- 



per limit other than our adopted one. However, we argue that 
our limits on the natal kick magnitude will not change signifi- 
cantly if we adopt a s malle r upper limit on the BH immediate 
progenitor mass (see {s!7. 11 1. 

Soon after the core collapse event, the BH companion fills 
its Roche lobe at periapsis. At this time, the BH companion 
has a mass and radius of 77.6lgn M Q and IOO9+377 Rq, re- 
spectively, and the binary orbital semi-major axis at periastron 
is 2186^02 ^0- At me en d of CE event, the binary consists of 
a 24.8-38.7M© BH and a 32.2-41.6M Q He star. These limits 
on component masses of the post-CE binary enclose the full 
range of our derived constraints, which are obtained from the 
modeling of the post-CE binary progenitor of IC 10 X-l (see 
® . Due to the intense stellar wind mass loss suffered by the 
He star, the binary orbital period continuously increases. In 
the meanwhile, a small fraction of the stellar- wind material is 
accreted onto the BH, resulting in X-ray emission. This pro- 
genitor eventually evolves to the observed XRB IC 1 X-l in 
< 0.36 Myr since the end of the CE phase (see jj7.2l >. 

We also find that the existence of IC 10 X-l gives very 
strong constraints on the physics of CE evolution. It is the 
very first time that an XRB appears to be capable of distin- 
guishing various treatments of CE evolution. A CE treat- 
ment is considered to be consistent with our derived post-CE 
properties of the BH-He star binary progenitor of IC 10 X- 
1, if it provides an CE effi ciency (ctrp) in the range of 0-1. 
We first follow the ori ginal IWebbinkl (119841) prescription and 
study the CE event, using models of the He star progenitor 
constructed with different initial mass es a nd v arious strength 
of convective overshooting (see 35.11 & 35.21 i. We find that 
acE is always > 1 for any models at any evolutionary age, 
because the envelope binding energy is always at least a fac- 
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tor of 3 larger than the maximum possible change of orbital 
energy during CE evolution. We then consider the formalism 
of hy per-c ritical accretion onto the BH during CE evolution 
(see {35. 3k In order to have olqe < 1, we find that the radius 
of the He star progenitor at the onset of the hyper-critical ac- 
cretion phase needs to be at least 1.3 times larger than that 
at the onset of CE phase. There cannot be any pre-CE bi- 
nary configuration satisfying this requirement. The last CE 
treatm ent considered in our a nalysis is the "enthalpy" formal- 
ism (Ivanova & Chaichenets 201 1). We find that it naturally 
provides chqe < 1 with realistic pre-CE binary configurations, 
because including the term P{m)/p(m) as an energy source 
lowers the binding energy of the envelope. Hence, the "en- 
thalpy" formalism seems to be the appropriate description on 
the energy budget of envelope ejection, at least for the mas- 
sive stars involved in the formation of IC 10 X-l. By adopting 
this CE treatment and considering all constraints on our de- 
rived evolutionary history of IC 10 X-l, we find acE to be in 
a range of 0.58-1 (at 95.4% of confidence). 

B ased on our derived evolutionary history of IC 1 X- 1 , the 
spin angular momentum of the BH immediate progenitor is 
likely to be low. This is because once it loses its H rich en- 
velope, the BH progenitor will suffer from the high mass loss 
rates of WR stars. This intense mass loss via stellar wind will 
take most angular momentum away from the BH progenitor 
and spin it down quickly. Tidal effects could have kept the the 
BH progenitor from spinning down. However, the pre-SN bi- 
nary orbit is relatively wide, with an orbital period of > 0.5 yr. 



This means that the tides exerted on the BH immediate pro- 
genitor by its companion star are expected to be weak. Even 
if the tidal interactions are much stronger than our expecta- 
tion, they can at best synchronize the spin of the BH imme- 
diate progenitor with the orbital frequency. Hence, under the 
assumption that spin angular momentum is conserved during 
the core collapse event, we argue that the natal spin of the BH 
in IC 10 X-l is likely to be small. Furthermore, since the BH 
has accreted negligible amount of mass from the stellar wind 
of its companion, it cannot be significantly spun up after its 
formation. Therefore, we expect the current spin of the BH in 
IC 10 X-l to be small as well. 

The BH spin of IC 10 X-l is expected to be measured in the 
coming years. If the spin of this BH turns out to be fairly 
high, it may imply that the BH was spun up by accreting 
significant amount of mass during CE evolution at a hyper- 
critical accretion rate. Another possible explanation could be 
that the BH was spun up duri ng the core collapse event (see 
iBlondin & Mezzacappa|[2007l) . 
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APPENDIX 
HYPER-CRITICAL ACCRETION ONTO BLACK HOLE DURING COMMON ENVELOPE EVOLUTION 

Let us denote the mass of the black hole (BH) by M B h, the mass of the BH companion and its core mass by M cgm and M core , 
and the binary semi-major axis by A. Accretion onto the BH will be initiated when the binary semi-major axis equals to the BH 
companion's radius. As the time interval between the onset of common envelope (CE) evolution and that of the accretion phase is 
relatively short, the masses of the BH and its companion are expected to be unchanged. Because of the very short circularization 
time scale due to the Roche lobe filling BH companion, the binary orbit is assumed to be circular at the onset of accretion. CE 
evolution and accreti on onto th e BH will end when the envelope of the BH companion is ejected to infinity. 

Following Belczvn ski et alJ (J2002), the rate of change in M B h and A with respect to M com due to accretion onto the BH is given 
by: 



[c d (M com +M BH )-M com ] 



dM B 
dM c . 



MnnMr, 



clA 



dM a 



■+M, 



BH, 



(Al) 



where Cd is the drag coefficient of the BH traveling in its companion's envelope. We adopt c c i to be 6 dShima et al.l 119851 : 
iBethe & Brown|[T998l) . In order to express these rates in ordinary parameters, we use the equation of energy balance during CE 
evolution. 
During the phase of accretion, the binding energy of the BH companion's envelope can be expressed as 



t^bind — Jej 



GM com (M com - M cor e) 

XA 



GM com (M com -M core - AM BH ) 
XA 



(A2) 



where f e j is the mass fraction of the ejected envelope as defined in Equation ( fTOb and AM B h is the amount of mass accreted onto 
the BH. Since the outcome of the CE event is known, which is given by the successful evolutionary sequences, we express AM B h 
using the known post-CE BH mass (M BH /) 

AM B H = M BH ,f-M BHl (A3) 

By substituting Equation dA3b into Equation dA2l ). we obtain 

GM com (M com - Moore ~ M BH j + M BH ) 



Eh, 



ind " 



XA 



(A4) 



The parameter A is a numerical factor scaling the binding energy of the BH companion's envelope during the phase of accretion, 
and is defined as 

LrlVL com.iv^ com,i "*core) , < c\ 

A = - • (A5) 

*- j bind,i I *~com,i 
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Here, M conv and R CO m,i are the pre-CE mass and radius the BH companion, respectively. We notice that R conh \ is different from 
the BH companion's radius at the onset of the accretion phase (R CO m.aa!)- The envelope binding energy Eund.i of the pre-CE BH 
companion is defined in Equation (|8). 

The envelope of the BH companion is ejected at the expense of the binary orbital energy. According to Equation (0, we can 
write 

a CE ■ -TT. = -TT. (A6) 

dM com dM c , 



L com ulri com 



Here, the rate of change in orbital energy with respect to M com is given by 



dM com 2A V ' dMcom A dM com 

From Equation (IA4b . we derive the rate of change in E b - m d with respect to M com as 



dEund G 



dM com XA 

Using Equation ( |A6b — dA8b , we can write 

2 \ dM BH 

1 + r \M com - 



dE orb G ( dM BH M BH M com dA 

W7 M com -JTZ ; 777 + M BH ( A7) 



(A8) 



dM BH M com (M core -M com +M bh . f -M BH ) dA 
M com — — + 'i — — + 2M com - M core - M bhJ + M hh 



acE^J dM, 



M a 



c<:>l n 



2 

M BH r (M core - M com + M BHJ - M BH ) 

a CE A 



dA 2 

M BH A2M com -M core -M BHJ +M BH ) (A9) 



dM com ' " a C £A 
Then, we can derive a system of two ordinary differential equations using Equation ( lAlb and dA9l ): 



dA A f l + ^) (A10) 



dM com M com \ acE^h 2 
dM BH 2M BH M, 



(All) 



dMcom a C E>^h\h2 
where 

hi = c d (M BH +M com )-Mcom (A12) 

h 2 = MbhMco "' (l + -2—\ +M BH + -L-(M com -M core -M BH . f +M BH ) (A13) 

«1 V a CE*J CtCE* 

Using the known post-CE parameters as initial conditions, we can obtai n the BH m ass (M B hj) and semi-major axis (A acc ) at 
the onset of the accretion phase by numerically integrating Equation ( lAlOt and ( IA1 lb fromM tw to M comi . In order to calculate 
the pre-CE semi-major axis (A,), we again use the energy balance equation acE^Eorb = &Ebmd- Here 

AE orb = \GM BH Mco m ,i (V t) (A14) 



and 



&E Mnd = ; -— ) (A15) 

^ \J<cc 



Since the radius of the BH companion equals to the binary semi-major axis at the onset of the accretion phase, we set 

Rcom,acc = A acc and obtain 

Z\lVl coln j iVi-core) I , **-acc 



Af- 



-t — v -■— i,uin,t ' - ^ Lwe/ i -, 



olce^Mbha \ Re 



(A16) 



As the binary orbit continues to shrink during CE evolution, A, > A acc . Also, the BH companion must completely fill its Roche 
lobe in the pre-CE binary orbit. Hence, A, needs to satisfy the constraint 

Rcom,i=Ai(l-edr L (A17) 

with a pre-CE orbital eccentricity (<?,) < 1. Here, r^ is the approximated Roche lobe radius (Eggleton 19831) . 
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